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Abstract 

It is shown, using a track-level simulation, that the use of the D~ —* K*°K~ 
decay channel for D~ reconstruction, in addition with the previously studied D~ — » 
mode, enables two fold gain in the ATLAS i^-mixing signal statistics through 
Bg — > D~a\ Bg -decay channel. A new modification of the amplitude fit method 
is suggested for the x s determination. Some general aspects of the S^-mixing phe- 
nomenon is illustrated by pictures of Casimir Malevich, Maurits Cornells Escher 
and Salvador Dali. 



1 Introduction 

The main features of the I? s -mixing studies were explained in previous notes (see 
also references cited therein), so we don't need to repeat them here. Instead we will try 
to give arguments that S s -mixing phenomenon is indeed worthy to be studied. 

The famous Russian painter Casimir Malevich said a long time ago: "The object in 
itself is meaningless ... the ideas of the conscious mind are worthless". We would like to 
choose his great painting "The black square", which is reproduced below, as a starting 
point of our introduction. 



But from this starting point it is possible to go to the very different directions de- 
pending from one's imagination. So let us imagine the following picture behind the black 
square |J : 

"A cat is penned up in a steel chamber, along with the following diabolical device: in 
a Geiger counter there is a tiny bit of radioactive substance, that perhaps in the course 
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of one hour one of the atoms decays, but also, with equal probability, perhaps none; if 
it happens, the counter tube discharges and through a relay releases a hammer which 
shatters a small flask of hydrocyanic acid. If one has left the entire system to itself for 
one hour, one would say that the cat still lives if meanwhile no atom has decayed. The 
state vector |\P > of the entire system would express this by having in it the living and 
the dead cat mixed or smeared out in equal parts." 

But this is of course nonsense, at least from cat's point of view! 

We have reminded Schrodinger's cat old story here in order to give an impression that 
although we all became familiar with particle mixing, because the superposition principle 
lies on a very background of quantum mechanics, this phenomenon is by no means obvious 
or trivial property of reality. 

But what is strange and queer at the macrophysics level can still appear as the most 
common thing at the microphysics level. It seems that even our existence is based on 
particle mixing as will be explained below. 

One of very important characteristics of elementary particle is its mass. We can get 
some insight about its origin from the following simple trick. The propagator of a massive 
fermion can be represented in such a way 

1 111111 
~~ — + ~ m - + -m-m- + • • • , 



p — m p p p p p p 
or graphically 

m m 
X X 

where a single line represents the propagator of a massless particle. So things look like 
as if the massless particle is propagating through some medium and the mass emerges 
as a result of friction or interaction with this environment. But what is the medium the 
particle interacts with? A (massless) fermionic particle can have the following interaction 
with some scalar field C int = gipipp: 

If now the self interactions of this scalar field are such that it doesn't disappear in 
a vacuum state and develops a nonzero vacuum expectation value <</?>, when it is 
convenient to expand ip =< p > +<//, where ip' corresponds to the physical scalar particles 
(excitations over the vacuum) and < ip > just gives the medium (the vacuum) where all 
of us are living. Now because of this decomposition of tp the fermion-scalar interaction 
splits into two parts: 

y v ^< # < y? > + y <p' 

The second diagram represents an emission of the real scalar quantum and the first 
one generates the fermion mass m = g < ip >. 
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But if, for example <i-quark can emit a scalar particle without changing its flavour, 
why can't it do this with changing the flavour? We know that the flavour is not always 
conserved, so the following interaction is not excluded: 



g \ cp - J<g <cp> + g) cp' 

d / df d 

But now the first term gives d — s mixing! As a result our initial d and s fields are no 
longer mass eigenstates (the states with definite mass), instead their time development in 
the rest frame is described by the matrix Schrodinger's equation (h = 1): 

_ ( m d m ds \ ( d\ _ 
\m ds m s )\s )- 

( cos 9 C sin 9 C 
y — sin 9 C cos 9 C 

where tan 29 c = m 2 "^ d and m d , fh s mass eigenvalues are defined from the equations m d = 
rhd cos 2 9 C + fh s sin 2 9 C , m s = rh d sin 2 9 C + fn s cos 2 9 C . It is obvious that the corresponding 
eigenvectors (the physical d and s quarks) are 






d = cos 9 C d — sin 9 C s 
s = sin 9 C d + cos 9 C s 



This particle mixing has one important observable consequence. If initially the weak 
transitions were possible only within the (u, d) or (c, s) pairs, now the intergeneration 
transitions u ^ s and c «-> d are also possible because, for example, the physical s- 
quark contains both d and s bare fields s = sin9 c d + cos# c s. Thus u — > s transition is 
proportional to sin# c - sine of the so called Cabibo angle. 

But we have three quark-lepton generations. So after the mixing the weak transitions 
are possible between any up and any down quarks. The amplitudes of these weak transi- 
tions are convenient to express as a 3 x 3 unitary matrix. This Kobayashi-Maskawa matrix 
is a generalization of the Cabibo angle and reveals a remarkable hierarchical structure M 



' V u d V us V u b 
V c d V cs V c b 
V V td V u V tb 



( 1-|A 2 
-A 

^ AA 3 (1 - p-irj) 




Here A = sin^ w 0.22 is a small quantity. So the intergeneration weak transitions are 
suppressed and this suppression is more strong for not neighboring generations. 

If rj 7^ 0, the Kobayashi-Maskawa matrix is complex and violates CP. It is commonly 
believed today that this CP- violation is an important ingredient in baryon-antibaryon 
asymmetry generation in the universe and so the source of our very existence. 

So far we were talking about particle mixing at quark level. But quarks are confined 
into hadrons and we can study quark-mixing only indirectly via hadron-mixing. 5-meson 
system is very promising in this respect: because of a large mass of the 6-quark we 
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can enjoy an asymptotic freedom advantage of QCD and calculate strong interaction 
corrections, unlike, for example, i^-meson system. 

In the Standard Model the Bd — Bd mixing originates from the following diagram (and 
from the second one there intermediate up-quark and W lines are interchanged) 

6 -^\AAA^ 



U; 



w 



w 



Un 



v, 



-d 



Ui stands for any up quark. So 



B d - mixing ~ ^ \\jl{m h rrij) 



where \ = V Ui bV*. d and I(mi,mj) represents the loop integral. This integral diverges 
quadratically. But this divergence is harmless because the unitarity of the Kobayashi- 
Maskawa matrix ensures its cancellation in the sum: the unitarity means \ = 0) 
therefore 

XiXjI(rrii, rrij) = ^ K^j [I( m h rrij) - 1(0, rrij) - I(m h 0) + 1(0, 0)] 



and these subtractions greatly improve the convergence. For example, in case of t-quark 
contribution, these subtractions lead to the replacement 



1 



1 



(k 2 



m 



2\2 



(k 2 



m 



2\2 



+ 



1 



k 2 (k 2 — m 2 ) k 2 k 2 k 2 k 2 (k 2 



m, 



From this expression it is also clear that in fact just t-quark contribution is dominant 
for i?-mixing, because of its extraordinary large mass. 

When ARGUS made his measurement of the ^-mixing ||, nobody thought that t- 
quark is so massive. So the result of this measurement appeared as a big surprise. We can 
even say that t-quark was discovered by ARGUS, because the large B d - mixing, observed 
by ARGUS, is very difficult to explain without the existence of the t-quark with mass 
> lOOGeU. 

We can infer from the above given diagram that even larger mixing is expected in 
-Bs-system: 



B s — mixing 



Bh 



mixing 



Is 



td 



1 



X 2 [(l - p) 2 + T] 2 



1-25 



We see also that the relative magnitude of B s and Bd mixings measures (1 — p) 2 +rf - one 
side of the notorious unitarity triangle 0. It is worthwhile to mention that this ratio is, 
to a great extent, free from hadronic uncertainties, which arise when we ask how quark 
and antiquark from the above given 5-mixing diagram really form 5-meson. 
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To summarize, the 5 s -mixing studies are interesting, because they reveal a very fun- 
damental underlying phenomenon - the generation of particle masses and mixing angles 
via the Higs mechanism, the least understood thing in the Standard Model. Because of 
heaviness of the 6-quark and asymptotic freedom of QCD, the theory gives very definite 
predictions about expected 5 s -mixing, hampered only from uncertainties due to our in- 
ability to solve QCD in the confinement region. But these uncertainties are also, to a 
certain extent, under control 0. The theoretical predictions involve such a fundamental 
property as the unitarity of the Kobayashi-Maskawa matrix (the existence of only three 
generations). Any deviation between the theory and experiment can lead to significant 
change of our present day picture of the elementary particle world (recall the ^-mixing 
story). The forthcoming ATLAS experiment sensitivity to the S s -mixing covers the Stan- 
dard Model prediction range || . So it will either give one more conformation of the theory 
or will open a window into a physics beyond the Standard Model. 



2 Black Square view on the B s -mixing 

We hope the above given considerations convinced the reader that the Black Square can 
hide a very reach content behind it. For example, reflecting about the Schrodinger's cat 
we can end with the following picture of -B^-oscillations (experimentally confirmed by 
ARGUS @) 



escherl.jpg here 
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As we have already mentioned, more rapid oscillations are expected in the B s system. 
So if man eagerly stares on the Black Square he will at last distinguish the picture of 
5 s -mixing: 



escher2.jpg here 



The most impressive people can remark even a symbolization of the CP- violation on 
this illustration: the nature for some reason makes an absolute difference between particle 
and antiparticle, not just the conventional difference between black and white. 

3 D~ — > K*°K~ decay channel for D s reconstruction 

To observe the I?!?-mixmg in a real experiment like ATLAS and extract the corresponding 
x s parameter, which characterizes the B s — B s oscillation frequency, you need to recon- 
struct B s meson and determine its decay vertex with great precision. Two of the B s decay 
channels were considered for this goal up to now: B° s —>■ D~n + |T[j and B° s — > D~di 0. 
For the second channel Dj — > 07r~, <fi — > K + K~ decay mode was used for the D~ recon- 
struction. It was mentioned in 0, || that other decay channels of D~ can be also used to 
increase signal statistics. In the present note we consider D~ — > K*°K~, K*° — K + tt~ 
decay mode as one of the possibilities: 
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I^tag 



As it is clear from the Table 1 below, this decay channel is quite promising if compared 
to the previously used Dj — > 07r~. 

Table I. 

Branching ratios and signal statistics for — > D~a>i(l2Q0). 



Parameter 


Value 


Comment 


L [cm~ 2 s _1 ] 


10 33 




t[s] 


10 7 




a(bb -> /iX) [fib] 


2.3 


p£ > 6 GeV/c 






\rf\ < 2.2 


N(bb -> /iX) 


2.3 x 10 10 






0.112 




Br(B° s -> Djaf) 


0.006 




5r(a^ -> p°7T + ) 


~ 0.5 




5r(p° -> 7T-7T+) 


~ 1 




Sr(D7 -> 07T-) 


0.036 






0.491 




Br(D~ -> K* ^) 


0.034 




firf^* -> X+7T-) 


~ 0.65 




Ar(ir + X-7r-7r + 7r-7r + ) 


136600 


D- -> 07T- 


A^(K+X-7r-7r + 7r-7r + ) 


170800 


£>7 -> x* u is:- 
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Event generation and reconstruction procedures are similar to ones considered in 0. 
All other general parameters, like impact parameter resolution for smearing and transverse 
momentum resolution, are also the same as in p| and can be found there. 

Contrary to the D~ — > 07r _ case, where D~ peak was clearly seen in the invariant 
mass distribution of the three properly charged particles, assuming that two of them are 
if-mesons and one is pion, now D~ peak is not seen (Fig. 1), perhaps because K*° is 
too wide as compared to <p. Thus the combinatorial background from the signal events 
alone is already able to hide the D s and the reader is left with the sentence || "if you 
should see the word "buffalo" written on a cage containing an elephant, don't believe 
your eyes". Although when D~ meson is reconstructed from its true decay products 
the resulting invariant mass resolution, shown on Fig. 2, is almost the same as for the 
D~ — > (p7T~ mode. The another picture on Fig. 2 shows K*°, reconstructed from its true 
decay products. 

The resolution in the _B s -decay proper time (Fig. 3a) o T rs 0.061ps is practically the 
same as for the D~ — > 07r~ mode. The corresponding £? s -decay length resolution in the 
transverse plane is ~ lOOyum and the relevant distribution is shown on Fig. 3b. 

We expect that signal to background situation when using Dj — > K*°K~ mode will 
be similar to what was considered in 0. Compare for example Fig. 4 from and from 
this work, which describes a possible background from the B® — ► D~af decay when D~ 
is reconstructed via D~ — > <pn^ or Dj — > K*°K~ modes respectively. 

The main reason which allowed a good signal to background separation in || was 
the fact that D~ and B® masses are shifted from the D~ and B® masses by about 100 
MeV. But this equally applies to the Dj — > K*°K~ case also, because our cut on the D~ 
invariant mass doesn't change very much. The only problem which can arise is a large 
K*° width (as compared to the 0-meson width) and therefore the cut on the invariant 
mass of K*° should be considerably loose. At present we are not aware of a background 
for which this circumstance will play a crucial role. 

Acceptance and analysis cuts are summarized in Table 2. As it is known f|Tl]], a 
second level trigger is necessary to reduce an event rate, which is still too high after the 
first level trigger (the tag-muon). For the B s — ► D~n + , — > 07r~ mode the problem 
was studied in |TIJ and it was shown that some loose cuts on the invariant masses of <fi 
and D s candidates can be used for this purpose. The resulting trigger efficiency appeared 
to be about 0.54. For the B s — > Djaf, D~ — > K*°K~ channel, discussed in this note, 
the second level trigger should be specially investigated, of course. However we expect 
that the similar mass cuts on D s and K*° candidates will work in this case also and a 
50% trigger efficiency should be a safe estimate. 
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Table 2. 

Number of signal events from B° s — > L>7 6^(1260) channel expected in ATLAS after 1 

year (10 7 s) of operation at 10 33 cm _2 s _1 . 



Parameter 


Value 


Comment 


N(K + K-1T-1T + 1T-1T + ) 


136600 


D- 


— > 07T 


N(K + K'7l~7l + 7l~7l + ) 

J. T 1 J. A. J. i. /I /I /I /I J 


170800 


Ds 


-> K*°K- 


Cuts : 

p T > 1 GeV/c 
\r}\ < 2.5 






A/"fX + X _ 7r _ 7r + 7r _ 7r + ) 


9015 (6 Q%) 


D- 




A^(ir + fs:-7r-7r + 7r-7r + ) 


9910 (5.8%) 


D- 




< 35° 
A0™ < 15° 

IM - M ol < 192 MeV/c 2 (±3a) 
M^tt - M a +| < 300 MeV/c 2 






Ai/i,,,, <r 1 fl° 

L -^H J KK ^ lu 
^&KK < 10° 

- M | < 20 MeVVc 2 
- M D - | < 15 Me\//c 2 




D- 


— > 07T 


Ay^ < 20° 

— A/f^ nl <r 8D MpV/r 2 
M„„ — M 1 <r 20 MpVIc 2 

| 1 VI J{ 1 VI — | \ 4J\J IV1 V J Ls 




D- 




N(K+K-7l-7l + 7l-7l + ) 


6830 (5.0%) 


D~ 


— > 07T~ 


N(K + K~Tl-Tl + Tl-Tl + ) 


6830 (4.0%) 


D- 




Dj vertex fit x 2 < 12.0 

rt~y vprfpv fit <^ 1 9 

fi? proper decay time > 0.4 r>s 
5° impact parameter < 55 \im 
B° s p T > 10.0 GeV/c 






N (K + K~ n~ n + n~ n + ) after cuts 


4100 (3.0%) 


D- 


— > 07T 


N (K + K~ n~ n + n~ n + ) after cuts 


4780 (2.8%) 


D- 




Lepton identification 
Hadron identification 
Trigger efficiency 
Mass cut ±2<j 


0.8 
(0.95) 6 
0.5 
0.95 




N (K + K~ n~ n + n~ n + ) reconstructed 


1240 (0.9%) 


D- 


— >■ 07T 


N {K + K~ n~ it + it~ it + ) reconstructed 


1330 (0.8%) 


D- 
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As we see, about 2570 reconstructed B® mesons are expected for 10 4 p6 _1 integrated 
luminosity from B° s — > D~af channel when both of D~ — > <f)7t~ and D~ — > K*°K~ 
modes are used for D~ reconstruction. To them we should add 3640 B® mesons from 
B® — > D~7i + channel M. So in total we expect 6210 reconstructed B® mesons per 10 4 pb~ 1 
integrated luminosity when all as yet considered decay modes are used. 



4 Peak amplifier 

To extract the oscillation frequency, from M.C. or experimental asymmetry distributions, 
the so called amplitude fit method is useful IT] . Here we describe some refinements of 
this method. 

In the amplitude fit method an asymmetry distribution A(t) is fitted with the cosine 
function Af it cos (x s t/r) in which x s is fixed and Af it is the only free parameter. Repeating 
the fit for different values of x s , we get Af it (x s ) distribution. This distribution is peaked 
at x s which corresponds to the true value of the oscillation frequency. 

The peak position in the Afi t (x s ) distribution can be found with the help of the 
recently suggested "quantum" peak finding algorithm |]I2"||. The idea of this algorithm is 



based on the property of small quantum balls to penetrate narrow enough obstacles. So 
if such a ball is placed on the edge of some potential wall it will find its way down to the 
potential wall bottom even if the potential wall is distorted by statistical fluctuations. 

Let us introduce instead of continuous x s some discrete parameter i, say through 
Ni = Afu(i/2). The transformation Afi t (x s ) — > u(x s ), which we call peak amplifier, is 
defined for selected discrete values of follows 



u(i/2) = Ui , u i+ i 



P; 



P 



-u. 



U; 



1 . 



(2) 



And P% t i±\ transition probabilities are determined by the initial iVj spectrum |12 



P M± i = AiJ2 ex P 

k=l 



N i±k - Ni 



i±k 



(3) 



1 condition, cr,- is a 



Ai normalization constant being defined from the Pi t i~\ + Pi,i+i 
standard deviation (error) of Ni as determined by the cosine fit. 

If now we apply this peak amplifier to the data after the amplitude fit we get the 
probability distributions shown on Fig. 5 (for x s = 30) and on Fig. 6 (for x s = 45). 
As we see, the peak amplifier enables a clear determination of x s from the amplitude fit 
spectrum. 



5 x s sensitivity range 

To estimate the ATLAS sensitivity range for the x s measurement, the analogous procedure 
was used as in 0. 

Amplitude fit is applied to the asymmetry distribution generated by Monte-Carlo 
program. The input parameters of this program, such as signal to background ratio, 
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Bg lifetime, proper-time resolution and dilution factors are the same as in 0, with the 
exception of the number of signal events, which was taken to be 6210. 

The amplitude fit spectrum is further transformed using the peak amplifier transfor- 
mation as described above. In the resulting u(x s ) spectrum the mean value of x s and its 
standard deviation is calculated considering probability density. The "experi- 

ment" is considered as successful if the measured x s value (mean value of x s according to 
the u(x s ) distribution) is within two standard deviations from the true x s -value defined 
in the Monte-Carlo program. 

For each x s point 1000 such "experiments" were generated and the fraction of the 
successful "experiments" was calculated. The highest value of x s , for which this fraction 
is above 95%, is considered as a sensitivity limit for the ATLAS experiment. This limit 
was found to be about x™ ax = 42.5. This is almost the same number as found in |§. In 
fact the peak amplifier method doesn't give a significant increase for the sensitivity limit, 
but it allows a more accurate x s determination, as is indicated by Fig. 6, because the 
probability peak is much more narrow. 

Fig. 7 shows the distribution of the x s values, found by the peak amplifier method, 
for 1000 "experiments", generated with the "true" x s = 42.5. 

6 "Where is the beginning of the end that comes at 
the end of the beginning?" 

So we are at the end of our investigation. Our main conclusions are: 

• D~ — > K*°K~ mode enables a two fold increase in the signal statistics for the 
B° s — > D~a\ decay channel. 

• the ATLAS experiment can reach a sensitivity limit for high 42.5 
with the certainty. 

• D~ — > K*°K~ mode can be used also for the B s — > D~tc + channel. If the same 
increase in signal statistics is assumed, the total number of reconstructed B s events 
can reach 10 4 per 10 4 pb~ 1 integrated luminosity. According to estimates from M 
this will mean a sensitivity limit for x s about 46. 

We began our story with the Black Square. Here is one more Black Square image 
which illustrates the continues progress in the i^-mixing studies. 
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escher4.jpg here 



This picture leads to a question about a meaning of the scientific progress, or even to 
a more general question |]13| "What profit hath a man of all his labour which he taketh 
under the sun?". We refrain to give any other comment about this things because flU 
"Wo alle Worte zu wenig waren, dort ist jedes Wort zu viel". 

Nevertheless we don't want to end with only black and white images of the i? s -mixing 
related stuff. With the very great imagination you can catch sight of the glorious full 
colour artist's view on the particle mixing behind the Black Square (note the role of 
vertexing in emergence of this dream): 
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tigers.jpg here 
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Farewell 

When we began this investigation our intention was to write a vivid and joyful story about 
5 s -mixing, some mixture of science and art. Only scientific framework appeared to us as 



too narrow to embrace the beauty of life, because |15j "All things are full of labour; man 
cannot utter it: the eye is not satisfied with seeing, nor the ear filled with hearing.". 

Unfortunately Sasha Bannikov suddenly died at the end of last year and we are forced 
to end this project without him. 

Farewell Sasha! Let this article be a small thing that remains after you in this world 
as your memory. 
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Figure 1: Invariant mass distributions of three charged particle combinations in signal 
events, assuming 2K + n (a) or 3n combination (b) . 
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Figure 2: Invariant mass distributions of reconstructed K*° and D s events. 
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Figure 3: Proper time (a) and transverse radius (b) resolutions for the reconstructed B\ 
decay vertex. 
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Figure 4: Six particle invariant mass distribution corresponding to the meson. Dashed 
line - expected upper limit for background from B° decay. 
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Figure 5: amplitude A(x s ) and probability u(x s ) for x s = 30. 
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Figure 6: amplitude A(x s ) and probability u(x s ) for x s = 45. 
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Figure 7: The distribution of the x s values, found by the peak amplifier method, for 1000 
"experiments", generated with the x s = 42.5. 
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